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Hydrodynamics can be formulated as the gradient expansion of conserved currents in terms 
of the fundamental fields describing the near-equilibrium fluid flow. In the relativistic case, 
the Navier-Stokes equations follow from the conservation of the stress-energy tensor to first 
order in derivatives. In this paper, we go beyond the presently understood second-order 
hydrodynamics and discuss the systematisation of obtaining the hydrodynamic expansion to 
an arbitrarily high order. As an example of the algorithm that we present, we fully classify 
the gradient expansion at third order for neutral fluids in four dimensions, thus finding the 
most general next-to-leading-order corrections to the relativistic Navier-Stokes equations 
in curved space-time. In doing so, we list 20 new transport coefficient candidates in the 
conformal and 68 in the non-conformal case. As we do not consider any constraints that 
could potentially arise from the local entropy current analysis, this is the maximal possible 
set of neutral third-order transport coefficients. To investigate the physical implications of 
these new transport coefficients, we obtain the third-order corrections to the linear dispersion 
relations that describe the propagation of diffusion and sound waves in relativistic fluids. We 
also compute the corrections to the scalar (spin-2) two-point correlation function of the third- 
order stress-energy tensor. Furthermore, as an example of a non-linear hydrodynamic flow, 
we calculate the third-order corrections to the energy density of a boost-invariant Bjorken 
flow. Finally, we apply our field theoretic results to the A/” = 4 supersymmetric Yang-Mills 
fluid at infinite’t Hooft coupling and infinite number of colours to find the values of five new 
linear combinations of the conformal transport coefficients. 
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I. INTRODUCTION 

Despite the fact that the history of research into phenomena involving the behaviour of fluids 
dates back millennia, hydrodynamics remains a subject of intense research both in mathematics 
and physics. This is a result of its extremely wide applicability across length scales: from the 
dynamics of quark-gluon plasma, to the evolution of the universe as a whole. In modern field 
theoretic language, the theory of hydrodynamics can be understood as the long-range, infrared 
effective theory that can be expressed in terms of the gradient expansion of the relevant fields. 
Recently, the formulation of hydrodynamics in the language of an effective field theory has been 
explored in several works, among them in [1-14]. 

In the non-relativistic limit, hydrodynamics is described by the Navier-Stokes equations, i.e. 
the equations of motion of energy and momentum transport, 

dop + V ■ {pv) = 0, (1) 

p{do + V ■V)v = —VP -|- -|- (C + ?//3) V (V • u). (2) 

Here, p denotes the sum of the fluid’s energy density e and pressure P. The two viscosities, the shear 
rj and the bulk encode the microscopic properties of the fluid. In view of the gradient expansion, 
these famous equations, that is the continuity equation (1) and the Navier-Stokes equation (2), 
can only sufficiently describe fluids at low energies. More precisely, the Navier-Stokes equations 
are sufficient to first order in the small parameter that controls the hydrodynamic approximation: 
kimfp, where k is the momentum scale and imfp the mean-free-path of the underlying microscopic 
processes. 

To describe fluids at higher energies, like the quark-gluon plasma being tested at the LHC and 
RHIC, we must work with the relativistic version of hydrodynamics. Furthermore, it is important 
that the theory permits for microscopic processes at higher momentum scales. Hence, it is natural to 
expect that higher-order corrections in kimfp will play an important role in the gradient expansion 
series, relevant for such high-energy fluid flows. This is the main motivation for this work. 

Second-order terms, which enter as the leading-order corrections to the Navier-Stokes equations, 
were first considered by Burnett [15]. In a more modern language, they were studied by Muller, 
Israel and Stewart [16-19], whose work was initially motivated by the well-known problem that 
first-order hydrodynamic solutions of the Navier-Stokes equations suffer from acausal propagation. 
They showed that the presence of second-order terms could cure these issues. In particular, the 
Miiller-Israel-Stewart theory is of phenomenological nature [20] where the variables are the energy 
density e, velocity and the viscous part of the momentum flow n“^, which are governed by the 
conservation laws = 0 along with the phenomenological equation = —H"^ — r]a°'^. 

The coefficient rn is known as the relaxation time, D represents the longitudinal derivative taken 
in the direction of the fluid flow, D = u°‘'S/a, and is the relativistic shear tensor, to be defined 
later. These equations are hyperbolic and encode non-instantaneous relaxation of momentum. 
Given that this is a phenomenological equation, one can naturally extend it beyond just linear 
order on the right-hand side. In fact, not only is this extension natural but it becomes necessary 
once conformal fluids are considered, in which case non-linear terms like • u) are present in 
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the stress-energy tensor. 

While second-order hydrodynamics has been important for stabilising hydrodynamic simula¬ 
tions, to our knowledge, no second-order transport coefficient has yet been measured. 

The Miiller-Israel-Stewart theory did not include all possible terms consistent with the sym¬ 
metries of hydrodynamics - a systematic procedure that is normally dictated within the effective 
field theory approach. Their work was extended in [20, 21], where the full second-order conformal 
hydrodynamic stress-energy tensor was constructed. In [20, 21], it was shown that such an exten¬ 
sion required five new transport coefficients, rn, k, Ai, A 2 and A 3 . In the non-conformal case, all 
possible structures were found by Romatschke in [22], resulting in fifteen transport coefficients. 

Construction of hydrodynamics normally asserts the existence of a local thermal equilibrium 
and an entropy current. Consistency with the second law of thermodynamics then requires global 
entropy production to be non-negative. It is well known that this condition imposes constraints on 
the values of the hydrodynamic transport coefficients: In first-order Navier-Stokes hydrodynamics, 
one finds that shear and bulk viscosities must be non-negative, r? > 0, C > 0. For a pedagogical 
discussion, see [23]. In higher-order hydrodynamics, entropy production is more subtle. Many 
recent works, among them [22, 23, 31, 33, 51-54], have investigated this issue. To find constraints 
on second-order transport coefficients, they imposed a much stricter condition - namely that each 
linearly-independent (tensorial) term in the divergence of the local entropy current must be non¬ 
negative, even those formally sub-leading in the gradient expansion. It was shown by Bhattacharyya 
[52] that under these condition, the number of independent, non-conformal second-order transport 
coefficients is reduced from fifteen to ten. The five conformal second-order transport coefficients 
remain un-constrained. While such a construction certainly suffices in ensuring non-negative en¬ 
tropy production, it is not completely clear whether it is necessary in all physical hydrodynamic 
systems. 

The question of entropy production in higher-order hydrodynamics become even more myste¬ 
rious if one takes into account the results of recent numerical holographic simulations combined 
with experimental observations of the dynamics of strongly coupled quark-gluon plasma. A robust 
and remarkable property that they display is that the hydrodynamic description of time-evolution 
becomes extremely accurate very shortly after the collision of two heavy ion nuclei.^ Typically, this 
hydrodynamisation time th is of the order of Ifm/c. In many cases, the hydrodynamisation time 
is shorter than the thermalisation time, which means that hydrodynamic evolution may become 
applicable even before local thermal equilibrium is established in the plasma. What does hydro- 
dynamic evolution in absence of a local thermal equilibrium imply for the existence of an entropy 
current? Furthermore, is there a notion of hydrodynamics at zero temperature? 

In this work, we will leave these important questions unanswered. As we will only study the un¬ 
constrained hydrodynamic gradient expansion, the new transport coefficients that we will find will 
correspond to the maximal possible set of transport coefficients in conformal and non-conformal 
fluids. 

While exactly conformal high-energy fluids are not known in nature, the beta function in asymp- 


^ For a recent review, see [24]. 
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totically free theories, such as QCD, becomes very small at high temperatures. In such regimes, 
conformal hydrodynamics can accurately approximate the behaviour of fluid-like states. Further¬ 
more, conformal hydrodynamics is extremely important for understanding the AdS/CFT corre¬ 
spondence, within which the connection between gravity duals and hydrodynamics was established 
in [25, 26]. As was clearly demonstrated by the duality, the IR limit of various large-quantum 
field theories behaves in accordance with the hydrodynamical gradient expansion. This is apparent 
for example from the dispersion relations of the extreme IR modes in the spectrum [27, 28], which 
can be expressed in a series expansion, 

Oh 

u = (3) 

71=0 

where an can depend on any of the m-th-order hydrodynamic transport coefficients, with m < n. 
The order of the hydrodynamic expansion Oh to which it is sensible to expand the series in Eq. 
(3) is presently unknown. The reason for this is that adding the terms above some order Oh 
may cause the sum to rapidly grow and spoil the series approximation to the mode’s physical 
behaviour. Leaving these issues aside for the moment, it is important to note that the most 
powerful feature of the gravitational dual of some field theory is that it encodes full information 
about the entire hydrodynamic series - the structure of the gradient expansion and the values of 
its microscopically determined transport coefficients. Beyond our desire to understand the general 
structure of relativistic hydrodynamics suitable for the description of high-energy fluid-like states 
of matter, AdS/CFT should therefore be in itself a very important motivation for studying higher- 
order hydrodynamics.^ 

The problem that we will study in this work is precisely how hydrodynamics can be systemati¬ 
cally generated at all orders in the gradient expansion, i.e. higher powers in kimfp- What is clear is 
that the complexity of such an expansion makes it rather forbidding for the higher-derivative terms 
to be analysed in the usual way. For this purpose, we will formulate a computational algorithm 
that can generate relativistic hydrodynamic gradient expansion at any order. To demonstrate 
its power, we will directly extend the works of [20, 22] and classify both conformal and non- 
conformal uncharged hydrodynamics at third order in four space-time dimensions. We will show 
that third-order hydrodynamic stress-energy tensor requires us to introduce 20 new tensorial struc- 
tures (and transport coefficients An with n = 1, 2, ..., 20) for conformal fluids and 68 tensors in 
the non-conformal case. Thus, we will find the most general next-to-leading-order corrections to 
the relativistic uncharged Navier-Stokes equations, with terms up to and including O (5^) in Eqs. 
(1) and (2). 

Next, we will use the AdS/CET correspondence to compute some of the conformal, third-order 
transport coefficients, in the supersymmetric = 4 Yang-Mills theory with an infinite number 
of colours Nc and at inhnite ’t Hooft coupling A. Using the linear-response dispersion relations for 
shear and sound modes, a two-point function of the stress-energy tensor and the non-linear boost- 
invariant Bjorken flow, we will derive the values for five linear combinations of these coefficients, 


^ For works on second-order hydrodynamics in AdS/CFT, see [20, 21, 29-36] and references therein. 
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From the point of view of phenomenological hydrodynamics, there is no clear way to define the 
basis of linearly-independent tensor structures at a given order in the gradient expansion. As a 
result, any linear combination of the transport coefficients can be declared as a set of independent 
transport coefficients. 

Beyond hydrodynamics to third order, recent works [37-40] have studied the linearised sub¬ 
sector of the hydrodynamic expansion to all orders. In this work, however, we will be interested in 
the full non-linear extension of the relativistic Navier-Stokes equations in all of the relevant fields. 

There are several other corrections to hydrodynamics that could in principle invalidate the 
gradient expansion at some order. Firstly, it has long been known that the correlation functions 
of conserved operators, such as (TabTcd), exhibit non-analytic behaviour beyond first-order hydro¬ 
dynamics. This effect is known as the long-time tails and was recently studied and reviewed in 
[23, 41, 42]. However, because these non-analyticities arise from loop corrections, they are sup¬ 
pressed by inverse factors of the number of colours in non-Abelian gauge theories (1/A"c) [41]. Thus, 
in large-A^c held theories, such as those most readily described by the gauge-gravity duality, long¬ 
time tails are sub-leading compared to the gradient expansion series. Furthermore, the interplay 
between the non-analyticities and the higher-order expansion in k has not yet been fully analysed 
and understood. 

Secondly, as already mentioned above, an important question regarding the hydrodynamic ex¬ 
pansion is the convergence of the series (cf. Eq. (3)). It is believed that the hydrodynamic series 
is an asymptotic one, analogous to the perturbative series in QFT. Recent works [43, 44] have 
studied the behaviour of the hydrodynamic expansion from this point of view and [43] showed that 
in the holographic Bjorken flow [45-48], as expected from a divergent asymptotic series with a 
zero radius of convergence, the series indeed breaks down at some order. However, the order at 
which the series may break down in general is unknown and therefore it is important to study the 
expansion at higher orders. Similarly, in QED, the results of the renormalised perturbation theory 
are expected to be divergent [49]. However, higher-order loop computations in QED have yielded 
some of the most successful high-energy theoretical predictions, consistent with experiments. 
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This paper is structured in the following way: In Sec. II, we will proceed with our analysis 
of the hydrodynamic gradient expansion by first constructing a systematic algorithm that can be 
used for its classification at any order. We will then use the outlined procedure to write down 
the stress-energy tensor for conformal and non-conformal relativistic fluids at third order in Sec. 
III. In Sec. IV A, we will look for the linear dispersion relations and find corrections to diffusion 
and the propagation of sound. Moreover, we will compute the scalar (spin-2) correlation function 
of the stress-energy tensor and the energy density corresponding to the Bjorken flow. These 
results will then be used in Sec. V to compute some of the new transport coefficients in the 
AA = 4 supersymmetric Yang-Mills theory by using the AdS/CFT correspondence. The appendices 
are used to prove a statement about the required tensorial ingredients in the classification of 
hydrodynamics, simplifications that occur in tensors that do not include any curvature tensors and 
to list all third-order scalars, vectors and two-tensors that enter the gradient expansion. 

II. SYSTEMATICS OF THE CONSTRUCTION 
A. Hydrodynamic variables and the generalised Navier-Stokes equations 

Let us begin our exploration of hydrodynamics by identifying the relevant hydrodynamic vari¬ 
ables and the equations of motion that govern the behaviour of fluids. In this work, we will focus 
on the uncharged relativistic fluids, although the systematics we outline here can easily be applied 
to the classification of the gradient expanded Noether currents that would correspond to additional 
conserved charges in the system. The steps of the construction that we employ will closely follow 
those presented in [23]. 

An uncharged fluid on a manifold with an arbitrary metric gab can be described by two near¬ 
equilibrium functions; the velocity field u“(x) and the temperature field T{x). In the presence of a 
Noether current J“, one also needs to promote the chemical potential to a space-time dependent, 
near-equilibrium field, /i(x). For present purposes, we only assume the existence of a conserved 
stress-energy tensor, which we write in the gradient expansion of the velocity and temperature 
fields, 

{u, T) + r(f) (On, or) +... + ( 0 ”n,..., 0 ”T) +... = TfJ; + (9) 

The conservation equation, 

= 0, (10) 

then provides dynamical equations (of motion) for n“ and T. It is worth noting that Eq. (10) 
gives us d equations in d space-time dimensions, which can be solved by the d — 1 independent 
components of n“ together with a single scalar field T (n“ is normalised to nan“ = — 1). However, 
for conciseness, we will only work in d = 4 dimensions in this paper, although our construction 
could be easily generalised to an arbitrary number of dimensions. 

Following [23], it is convenient to decompose into 

= Su^v!’ + r (11) 
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where E and V are scalars, a transverse vector and a transverse, symmetric and traceless 
(TST) tensor. Here we have introduced the projector, 

(12) 

Each one of the operators £, V, q°‘ and is then gradient-expanded in rt“, T and gab, to facilitate 
the expansion in Eq. (9).^ 

In order to construct a physically sensible theory of hydrodynamics, it is not sufficient to only 
hnd all possible tensor structures for E, V, q°‘ and The first issue arises from the fact that we 
have not microscopically specified the meaning of the velocity and temperature fields. Because we 
are working with the gradient expansion, it should therefore be possible to transform the fields by 
adding to them terms sub-leading in the gradient expansion while keeping the physical predictions 
of hydrodynamics invariant, 

+ {dT, du) + {d^T, dTdu, d\) + ..., (13) 

T^T + {dT, du) + f^ 2 ),T {d^T, dTdu, d\) + ... . (14) 

The hydrodynamical description of a system should thus remain invariant under such frame re¬ 

definitions. To study uncharged fluids, it is most convenient to work in the Landau frame with 

= 0, (15) 

which requires us to set all = 0 and = 0, forn > 1. The zeroth-order = 0 because there 
are no transverse vectors without derivatives among our hydrodynamics variables. The remaining 
term in the expansion of E, i.e. E = T(o)) is then identihed as the energy density of the fluid, 

^ = ^(0) = After such an adjustment of the arbitrary functions in Eqs. (13) and (14), we are 

left with 

(16) 

where = P is the pressure of the fluid. Hence, what remains to be gradient-expanded in the 
classification of uncharged fluids, are the scalar V and the transverse, symmetric and traceless 
There is an additional source of restrictions that we must impose on the structure of the hy¬ 
drodynamic tensors in order to obtain the set of all “physically” independent structures: when 
working at the n-th order of the gradient expansion, we can always use the equations of motion 
(the conservation equation (10)) to eliminate some terms at the expense of introducing corrections 
at all sub-leading orders in the gradient expansion. More precisely, Eq. (10) gives 

= —Ds — {s + P)^ ■ u + higher derivatives = 0 , (17) 

A^VcT'^^ = {e + P) Du°^ + Vj_“P -|- higher derivatives = 0 , (18) 

where we have dehned the longitudinal and the transverse derivatives as 

D = u'^Va , , (19) 


® Throughout this work, we will be using the (—, -I-, -t) Lorentzian signature for the metric tensor. 
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so that Va = V_La — UaD. The most important property of the transverse derivatives is that when 
acting on an arbitrary tensor, 


u^'^±cVbib2... = 0 - ( 20 ) 

Eqs. (17) and (18) are the generalisations of the relativistic Navier-Stokes eqnations for fluids 
with a gradient expanded stress-energy tensor to an arbitrary order. Following [22], it is convenient 
to use the thermodynamic relation e + P = sT, where s is the entropy density, to re-write Eqs. 
(17) and (18) as 


In s = —V • u + higher derivatives , (21) 

Du'^ = —Cg V_l“ In s -|- higher derivatives . (22) 

The constant Cs is the speed of sound in the fluid. Note also that V • n = V_l • tt. As a result of 
the form of Eqs. (21) and (22), we will use the scalar function Ins instead of the temperature field 
T to construct the gradient expansion. 

The key observation that follows from Eqs. (21) and (22) is that it is most convenient to only 
work with the transverse derivatives of and Ins, and the Riemann tensor Rabcd- This is because 
all longitudinal derivatives of and Ins can always be written as purely transverse derivatives 
and commutators of the covariant derivatives, which can then be expressed in terms of various 
combinations of the Riemann tensor, its covariant derivatives and transverse derivatives of and 
Ins. The proof of this statement is presented in Appendix A. 

B. Constructing the gradient expansion 

We are now ready to begin classifying the hydrodynamic expansion. By choosing the frame for T, 
or the entropy density field s, and defining V(o), we have already fixed zeroth-order hydrodynamics, 
which describes the behaviour of ideal, non-dissipative fluids. There exist no transverse vectors, 
neither TST tensors that could be constructed out of only or gab- 

We can thus turn our attention to the hydrodynamic tensors with non-trivial derivative struc¬ 
tures. To construct the hydrodynamic expansion at higher orders, it is most efficient to carefully 
select the relevant ingredients that will build the irreducible tensor structures, i.e. those that can¬ 
not be eliminated by the equations of motion (21) and (22). As discussed at the end of Sec. IIA 
and Appendix A, all the longitudinal components of Vattfe and Va In s can be eliminated. There¬ 
fore, it suffices to use transverse derivatives of and Ins, as well as the Riemann tensor Rabcd and 
its covariant derivatives. 

Let us begin the classification of hydrodynamics by considering first-order tensors. The only 
possible one-derivative ingredients that can be used to build £, P, and are V±a'Ub and 
V_Lalns, as there are no one-derivative tensorial structures that would include the metric tensor 
Qab- We can therefore write the first-order set of tensorial ingredients as 


= {V_La«6, V_La In s} . 


(23) 
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Similarly, we define to be the set of all n-derivative objects relevant at the re-th order of the 
gradient expansion. The two and one-index tensors in Eq. (23) must then be contracted in all 
possible ways with the zeroth-order structures, i.e. Ua and Qah-, 

gab] . (24) 


Imagine that we wish to form all two-tensors with single entries from and an arbitrary 
number of the X^^^ structures. It is easiest to count such combinations by using covariant and 
contravariant indices on the two different sets, as it would be redundant to add additional strings 
of X^^) and then contract them amongst themselves. The total number of free indices from the 
tensors is ai -|- 2 a 2 , which should be directly contracted with either a 2 - or a 1 -index object from 
Hence, ai -|- 202 — 2 or oi -|- 202 — 1 can equal either +2 or — 2 . 


r(n) 


the 


To generalise the preceding discussion to tensors of any rank, let us now denote by 
number of un-contracted indices in the m-th entry of . All possible contractions of X^^'> with 
x(i) can thus be found by solving 


^o„ [Xfj - Mb 


= ±fV, 


(25) 


n=l 


for each m, such that o„ £ Z+U{0}. The three different values of = {0,1, 2} correspond to scalars 
(terms in £ and V), vectors (terms in q“) and tensors (terms in t“^), respectively. Throughout this 
work, we will use the symbol to denote the set of all possible tensors with n derivatives that 
can be used to form the tensors of rank N. 

Using Eq. (25), we find that in first-order hydrodynamics. 


m = 1 : oi -|- 202 — 2 = ±A^, (26) 

m = 2 : oi -|- 202 ~ 1 = (27) 

Eor scalars with = 0, the possible combinations are 

m = l: {( 01 , 02 )} = {(2,0), (0,1)}, (28) 

m = 2: {( 01 , 02 )} = {(1,0)}, (29) 

which give us the possible tensors, 

= {UaUh'^LcUd, OaV^follls} . (30) 


For a vector with N = 1, we need to solve for the right-hand-side of (25) with -XN = ±1, giving 
us a set of equations 


m = 1 : oi -|- 2 o 2 = 3 V oi -|- 202 = 1, (31) 

m = 2 : oi -|- 2 o 2 = 2 V oi -|- 202 = 0, (32) 

which have solutions 

m = l: {( 01 , 02 )} = {(3,0), (1,1), (1,0)}, (33) 

m = 2: {( 01 , 02 )} = {(2,0), (0,1), (0,0)}. (34) 
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The possible tensors from which the one-derivative vectors can be constructed are thus 

= {UaUbUcV ±dUe, UaQbc'^ UaS/^b^c, UaUb^ ±c Qab^ ±c In S, V_La In s} . (35) 


Similarly, the two-tensors in can be constructed by solving 


m = 1 ; ai -|- 2a2 = 4 V oi -|- 2a2 = 0, (36) 

m = 2 : ai -|- 2a2 = 3 V oi -|- 2a2 = —1, (37) 

which gives 

m = l: {(ai, 02)} = {(4,0), (2,1), (0,2), (0,0)}, (38) 

m = 2: 1(01,02)} = {(3,0), (1,1)}. (39) 

Hence, 

7^2 — {uaUbUcUd^ ^ UaUbQcd^A-e^f^ 9ab9c(^ A-e^fi ^-La^6) 

UaOftOcVxd In S, Oa5fecV_Lrf Ins}. (40) 


All the tensors in the sets and must now be contracted in all possible ways with 

Qab to form scalars, vectors and two-tensors, respectively. The next step in the construction is the 
elimination of duplicates and simplifications that use the properties of the participating tensors. 
In particular, we can use the fact that 


UaU 


(41) 


hence OaVj_;,o“ = 0, and so on. Furthermore, the property of the transverse derivatives that was 
stated in Eq. (20) can also be used to vastly simplify the resulting expressions. In fact, it should 
be clear that the combination of Eqs. (20) and (41) makes any contraction between in and 
the entries of vanish. However, this will not be the case for the curvature tensors in 

All scalars that are made from J'^'^ can immediately be used in S and V. On the other hand, 
the vector and tensor require additional treatment. Since must be transverse, we may 
only keep the vectors that vanish after they are contracted with Ua- To construct all TST tensors 
relevant for it is easiest to form manifestly TST structures out of all two-tensors by applying 
the operation 


1 




/^ab^cd^ 


cd • 


(42) 


All linearly independent tensors that survive this operation can then be considered as terms in 
Although the above procedure yields the correct results, it is computationally more efficient to 
construct the tensors that do not include any Riemann tensors separately. As already noted below 
Eq. (41), this is because our choice of the tensorial ingredients (Appendix A) vastly simplifies the 
possible contractions between X^^^ and the tensors excluding various curvature tensors in X*'”’^ with 
n > 1 (curvature tensors can have covariant derivatives acting on them). In Appendix B, we prove 
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the claim that when constructing scalars, transverse vectors and TST two-tensors with only 
and Ins at any order of the gradient expansion, we can only use the gab component of i.e. set 
fli = 0 in Eq. (25). For tensors containing Rated, we have to allow for any ui. 

Once all the linearly independent tensor structures are found, we associate hydrodynamic trans¬ 
port coefficients with each one of the tensors and write them as a sum in the stress-energy tensor. 
The transport coefficients cannot be determined phenomenologically, but must be calculated from 
the underlying microscopic theory.^ To denote the transport coefficients at all orders, we will use 
the lower case Greek letter “upsilon” with three indices. 


(n,N) 

i 1 


(43) 


where n is the order in the gradient expansion and N determines whether a particular transport 
coefficients comes from the scalar (V), vector (q“) or two-tensor (t“^) sector of r“^, denoted by 
= 0, = 1 and N = 2, respectively. In case one wanted to work in a frame with £ 7 ^ , some 

additional notation should be introduced to distinguish between the two cases. Finally, i 

will count the number of transport coefficients in each of the sectors, at each order of the gradient 
expansion. 

Since we are working in the Landau frame and in the frame with £ = T(o) = e, we only need to 
find 'P(i) and at hrst order. What we find from the above considerations is 

^(1) = -CV_L • U, 

fab _ „^ab 
^(1) — —, 

where = C and = g are the two hrst-order transport coefficients 

viscosity, respectively. The only hrst-order TST tensor is dehned as 

= (46) 


(44) 

(45) 

: the bulk and shear 


For future use, we also introduce the most general decomposition, 

^A-a^b — '2'^ab ^ab T ) (^^) 

where the vorticity is 

^ab ~ 2 (^-La'^6 V±b^a) • (^8) 

At second order in the gradient expansion, the non-linear tensorial ingredients are 

XG) ^ {V_LaV_Lfettc, V±aUb V_Lc«d, ±aUb V_Lc In S, V_LaV_L6 hi S, V_La In S V_Lb In s. Rated} , (49) 

where includes the Riemann tensor that involves two-derivative metric structures. The full 
structure of second-order hydrodynamics was hrst found by Romatschke in [22] and is given by a 


^ Such calculations can be rather involved and use a combination of linear response theory and kinetic theory 
techniques or lattice gauge computations. See for example [31, 50]. 
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sum of fifteen linearly independent tensors with fifteen transport coefficients, 


1 P(2 ) = (V • u) + + 6 (V • uf + + e4Vxa In In s 


ab{ 


+ ^bR + ^eUaUbR^^, 


(50) 


t$) = 


(V • u) 

O 


+ K 


+ \riT^a'^^ (V • u) 

O 


+ 2K*UcR^^‘^^^^Ud + + A4 Vx<“ In Ins. (51) 


For historical reasons, in the Miiller-Israel-Stewart theory [16-19], the transport coefficient r/rn 
is written as a product of the shear viscosity r/ and the relaxation time rn- Furthermore, the 
physically motivated combinations of the tensors (shear tensor, vorticity) that appear in Eqs. (50) 
and (51) are not necessarily those that would directly follow from our systematics, but can easily 
be written in terms of our tensors by simple linear combinations. What is invariant, regardless of 
the choice of the linearly independent tensors, is their number. 

For calculations that do not employ the Landau frame, it is useful to also write down the 
transverse vector q°‘. At first order, the only term that enters into the expansion is 

<i) =^^^'^Vx“lns, (52) 

while at second order, the five independent transverse vectors give 

<2) + uf ^V'’“Vxblns 

+ Ins + (53) 


Beyond the possibility to work in a different frame, the knowledge of the vectors is also required 
for building the gradient expansion of the entropy current 3°“ [22, 23, 51, 52]. However, in that 
case, we would not impose that 5“ be transverse. In this paper, we do not consider the entropy 
current analysis nor do we study potential constraints that may arise from demanding positive local 
entropy production [22, 23, 31, 33, 51-54]. As discussed in the introduction, it is well known that 
in first-order hydrodynamics this analysis results in the condition that shear and bulk viscosities 
must be non-negative [23]. In second-order hydrodynamics, [52] argued that the entropy current 
considerations give hve constraints, resulting in the reduction of the 15 second-order transport 
coefficients to 10 independent ones. 

In a different frame for the temperature field T, which would result in ^ ^ 7^ 0) 

the structure of the gradient expansion for would be exactly the same as for 'P(n)- In such 
cases, we would therefore only need to associate a new transport coefficient with each term in , 
independent from those in V{n). 

Finally, we note that by following the above procedure, we could also include the fluctuating 
chemical potential //(x), which would then allow us to construct the gradient expansion of the 
conserved Noether current, such as for example the baryon number current. We will not explore 
this option in this paper. 
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C. Conformal hydrodynamics 


A very important sub-class of fluids are the conformal fluids, which have played a central role 
in applications of the gauge-gravity duality to hydrodynamics [25, 26]. Fluids with conformal 
symmetry are characterised by a vanishing trace of their stress energy tensor, 


n = 0, (54) 

and the Weyl-covariance of T“^, i.e. each term in the stress-energy tensor needs to transform 
homogeneously under the non-linear Weyl transformations [20]. 

In quantum field theory, the perturbations of the metric tensor, gab —)• gab + hab, can be used to 
source the stress-energy tensor correlation functions. To first order in the source hab, the generating 
functional includes the coupling 

J d^X^T^'^hab. (55) 

Under the Weyl transformation, the metric tensor transforms with conformal weight of = —2, 

9ab e~‘^‘^gab ■ (56) 


Eqs. (55) and(56) then imply that in order for the stress-energy tensor to scale homogeneously, it 
must have the conformal weight Arpab = 6 in four dimensions. 

Tab Tab. (- 57 ) 

In general, a tensor •••“’" with the conformal weight = [A] + m — n that transforms homo¬ 
geneously under the Weyl transformation, transforms as 

A(l\...CLm _ A(l\...CLm 

^bl...brb - ^ ^bl...bn • 

In the expression for the conformal weight, [A] denotes the mass dimension of the tensor operator 
and m and n count the number of contravariant and covariant indices, respectively. Hence, Ua, Ins 
and the transport coefficients transform as 



^ (59) 

Ins —)• In (e^'^s) , (60) 

^{n,N) (61) 

For a detailed discussion of the Weyl transformations in hydrodynamics, we refer the reader to 

[ 20 ]. 

The procedure for finding the most general conformal hydrodynamic stress-energy tensor is then 
the following; since we are working in the Landau frame with only 7 ^ Oj f^ie traceless condition 
(54) implies that e = 3P. Furthermore, it implies that all with n > 1 must vanish. We 
must then use the full series of tensors inside the non-conformal expression for the traceless 
discussed in Section IIB, and hnd how each one of them transforms under the Weyl transformation 
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(58). Since (57) must be ensured, those tensors that transform homogeneously can be immedi¬ 
ately used in the conformal Those that do not, must be combined by linear relations into 
the maximal possible set of homogeneously transforming tensors. Each of the remaining linearly 
independent homogeneously transforming TST two-tensors can then be assigned an independent 
transport coefficient and the full series constitutes the stress-energy tensor of a conformal fluid. 

Up to second order and in the Landau frame, the above procedure gives us the conformal 
stress-energy tensor [20, 21], 


n ab j.ab , j.ab 

= ^( 1 ) + ^( 2 ) = -w + vm 


(V • u) 
o 

+ Aia<“ + A2 (t<“ 12^)^ + A3l2<“ 


-I- K 




(62) 


where 11“^ was dehned in Eq. (9). What remains from the list of all non-conformal transport 
coefficients is one first-order transport coefficient r] (shear viscosity) and five second-order transport 
coefficients rn, At, Ai, A2 and A3. 

Einally, we should note that we can completely ignore the effects of the conformal (Weyl) 
anomaly on the hydrodynamic transport in this work, as we are only considering third-order 
hydrodynamics in d = 4 dimensions. It is well known that the Weyl anomaly in four dimensions is 
given by 

/rna \ _ ^ ( p nabcd a p nab , p2'\ , ^ ( p nabcd tip nab , ^ p2 1 fa‘ 3 \ 

a) IGtt^ \^^o,bcd^ ^ j ' [^abcd^ ^^ab^ ' 3^ / 

~ {dgab)'^ ■ (64) 


Hence, this Weyl trace anomaly could only affect fourth-order hydrodynamics [20]. In that case, 
the use of only in the expansion of would be insufficient, as is manifestly traceless. In 
fact, Eq. (63) could simply be incorporated into in the Landau frame, by writing 


T> — ^ ^ f? nabcd . 2o C C So 2 , 

The three new transport coefficients. 


,(4.0) 

c — a 

’1 

487r2 ’ 

,(4.0) 

2a — c 

2 

247r2 

,(4,0) 

c — 3a 

3 

1447r2 


(65) 


( 66 ) 

(67) 

( 68 ) 


are thus completely determined by the central charges a and c of the underlying conformal field 
theory. We defer a more detailed study of fourth-order hydrodynamics to future work. 


III. THIRD-ORDER HYDRODYNAMICS 

We are now ready to construct the hydrodynamic stress-energy tensor at third order in the 
gradient expansion. To use the systematics described in Sec. II, we first need the set of third-order 
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tensorial ingredients, 

X(3) ^ {V±aV±bV±cUd,V_LaV±bUcV±dUe,V_LaUbV±cUdV±eUf, 

V±aV±brrc V_Ldlns, V_LaUbVxcV_Lrflns, VxaVxfeVxclns, 

V±aV±b In S V_Lc In S,V±a In S V_Lfe In S V_Lc In s, V_LaWb Rcdef, 

^^ R-bcde^^aRbcde} • ( 6 ^) 

We can then show that the third-order non-conformal stress-energy tensor, in the Landau frame, 
takes the form 


r“3')=P(3)A“' + tg), (70) 

where P(3) is a sum of 23 terms and a sum of 45 TST two-tensors, each with a new transport 
coefficient, 


23 45 

i=l i=l 

The scalars and the two-tensors are listed in Appendices D and C, respectively. In total, these give 
68 new transport coefficients. If we were to work in a different frame, it is also useful to list all 
possible hydrodynamic vectors that may enter the expansion of . The 28 transverse third-order 
vectors are listed in Appendix E. 

Out of the two-tensors, we can now construct conformal third-order hydrodynamics, as described 
in Sec. IIC. The expansion of the Weyl-covariant conformal stress-energy tensor takes the form 

20 

rp, = ET’o“‘, (72) 

i=l 

in which there are 20 new Weyl-covariant tensors that require us to introduce 20 new transport 
coefficients. In terms of 7)“^, the conformal tensors are 


^ 27 ii 7 i 2 27 i 3 27 i 4 27 i 6 27 i 7 Tis Tig T20 T2S Tao STiu 

'-^1 ~ ' 1 H-5-5—I-5-5-5-1--5-5-7-4/27-1—(7-7;—I-7;— 

« 3 ooo« 39 « 3 ob 1 Z Z 


+ 742 + 


3743 , 3744 


745 

2 ' 


(72-r2 + —+ —- —-Tie 


2Ti7 2 T 8 4T9 To j- Ts 3T9 3 To Ts 

~9 9 9 ^ ^ 2 ~Y 


Ts Ti 

2 2 


■ T2 — 


T3 , 5T4 


Ts 

2 ’ 


n -T 4 - r 

(P 3 _r 3 + —+ —-T 4 - 3 


2 2 

Ts 2T7 2 T 8 4T9 To T 8 To 5 To Ts Ts 

9 9 6 T “2 ^ “2 ^ 


(73) 


(74) 


9 


. Ti ^ Ts 3 T 4 
(P4_r4 + —+ —- —-Te 


Ta5 


2T7 2 T 8 4T9 To Ts ^ 


Ts Ts , T- 

~2-^ ~ 2 


(75) 

STi 


STs 


T4 

2 


Ts 

2 ’ 


^ 2 Ti Ts ^ Ts 2T7 2 T 8 4T9 To Ts 

( 7 s-Ts + — + — -T 4 - — + -g-g-g-— 


-Ts 


Ts Ts 


30 


To — 


(76) 

5Ti 


9 


9 


9 


2 
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Oe 


O 7 

Os 

O9 

Oio 

On 

O 12 

Ol 3 

Ol 4 

O 15 

O 16 

Oi7 

O18 

Ol 9 

O20 


=r6 


+ 


3r43 

2 

+ rii- 

3744 , 


= T7-i%- 

= Ts-%- 

T 20 

-rio + -^ 

= 721 — 372 


744 
4ri5 

3 

745 

7i7 _ 
3 

hi 

3 


745 

, 2ri7 


2 X 19 


720 


728 729 


3730 3733 


T- 720 , 3728 OT- 

■ /19- Y H-^-J /29 


Zl£ 

3 


3r27 - 


3729 3730 


3728 3740 , 9741 


3742 


3740 _ 9741 _ 9743 
2 2 2 

— 3743 + 3744 , 

9743 , 9744 3745 


28; 


= 722 — 3727 — 728 H—^ - 

= 723 ~ 728 ~ 729 j 

_ -r T ^h29 3 730 

/24 / 2 s “T 2 2 ^ 

= 725 ~ 728j 
= 726 ~ 2727 + Ts 


3r3 


30 


301 


=r3. + f 


T- 735 Tse 


^ + r45, 


2 X 35 


4736 2740 


4T45 


-%2 + 2%,- 3 3,3 

= 737 — 2740 + 6741 + 6742 + 6743 + 6744 + 27451 

= Tss — Tio + 2741 + 742 + 2744 , 

= Tsg — Tio + 2741 + Ta 2 + 2743. 


^ + 2r40-^ 


9744 

2 


3745 

2 


(77) 

9X43 

2 

(78) 

(79) 

(80) 

(81) 

(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 

(90) 

(91) 

(92) 


It should be understood that all O and T tensors carry two indices, a and b. Clearly, the choice of 
O is not unique and we could have chosen a different linear combination. What is important is the 
number of linearly independent combinations that gives us the number of corresponding transport 
coefficients. 

Furthermore, it is interesting to note that in flat space, only the first 15 Weyl-covariant tensors, 
Qab _ remain non-zero. This is consistent with the results from kinetic theory, within the 
context of which [55-57] found 14 transport coefficients in third-order hydrodynamics. Similarly, in 
second-order hydrodynamics, kinetic theory gave only 3 instead of 4 conformal transport coefficients 
in flat space. There, A 3 = 0, which is the transport coefficient corresponding to the term consisting 
only of vorticity in the gradient expansion, i.e. It is therefore plausible that the tensor 

missing from the kinetic theory analysis is the one consisting solely of vorticity tensors, i.e. ~ 03.5 

Finally, we stress that the number of new independent transport coefficients should be under¬ 
stood as an upper bound on their actual “physical” number (a maximal possible set), since there 
could be further constraints coming from imposing non-negativity of the entropy current diver¬ 
gence [ 22 , 23, 31, 33, 51-54]. As we have not yet performed the local entropy current analysis in 
third-order hydrodynamics, we cannot establish whether the non-negativity of entropy production 


® We thank Amaresh Jaiswal for bringing the results of kinetic theory to our attention. 
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imposes any constraints on the 68 and 20 un-constrained transport coefficients in non-conformal 
and conformal hydrodynamics, respectively. It is interesting to note that while such an analysis 
restricts the number of non-conformal coefficients, at least in second-order hydrodynamics, the 
number of conformal coefficients remains the same [52]. 

IV. PROPERTIES OF LINEAR AND NON-LINEAR THIRD-ORDER 
HYDRODYNAMIC TRANSPORT 

In this section, we analyse the effects of third-order hydrodynamics on linear transport as well as 
an example of non-linear hydrodynamics. First, we compute the dispersion relations for shear and 
sound modes to third order in the hydrodynamic expansion. Then, we compute the scalar channel 
stress-energy tensor two-point function. Since the dispersion relations and two-point functions 
follow from linear perturbations of in the hydrodynamic variables, only the stress-energy tensor 
terms linear in u“. Ins and gab are relevant. Lastly, as an example of a hydrodynamic process that 
is sensitive to non-linear terms, we compute the energy density of the boost-invariant Bjorken flow 
to third order. 


A. Shear and sound dispersion relations 


To compute the shear and sound dispersion relations in flat space, we need the linear part of the 
non-conformal third-order stress-energy tensor in Eq. (70) that includes derivatives of the velocity 
and entropy density fields. This is given by the following series of terms: 

Fay- = + ET"T“‘- (93) 

i=l i=l 

It is then easy to find the relevant dispersion relations by following e.g. [22, 34]. What we need to 
do is perturb and In s around the equilibrium values u°‘ = (1, 0,0, 0) and In s = const., and solve 
the generalised relativistic Navier-Stokes equations (21) and (22), that result from VaT“^ = 0, for 
uj{k) in the small k expansion. 

In flat space, V_La = UaU^db + da, which evaluated on the background, = (1,0,0,0), gives 
= Uado + da- Hence, to zeroth order in linear perturbations, Vj_o = 0 and Vj_j = di- For the 
purposes of linearised hydrodynamics, therefore 

d±a = (0, di) -|- non-linear , (94) 


which makes it clear that the linearised Navier-Stokes equations have no additional time derivatives 
acting on the hydrodynamic variables that come from higher-order hydrodynamics. Using the 
notation of Eq. (94), the linearised stress-energy tensor (93) in flat space becomes 


pab 

■ (3),/m 


^(3,0) ^ ^(3,0) ^ ^(3,0)^ 7]‘^^d±dd±'^d±aU'^ + (uf+ uf 


+ 


C 


(95) 
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In the two different channels, i.e. the transverse (shear) and the longitudinal (sound), we need 
to turn on the following fluctuations: 


shear: 

sound: 


5uy = 


_ —iLut-\-ikx 


<51ns = e-*‘^™5^,fc, (96) 

(51ns = e-“^^5^,fc. (97) 


The two dispersion relations in third-order hydrodynamics then follow from solving the equations 
of motion (21) and (22) with the linearised stress-energy tensor (95). We find that the dispersion 
relations, which are expressed in the form of Eq. (3), involve nine new third-order transport 
coefficients. 


shear: 


uj = —i 


V 


e + P 


— i 


1 V 


'2 + 2 


( 3 , 2 ) ,( 3 , 2 ) 


+ V. 


{e + PY 


e + P 


sound: w = +Csk — iVk^ ^ 


1 

2c. 


- 2c? 


2 ??rn 1 Ct 


3e + P 2e + P 


( 3 . 2 ) 

^ k^ + O {k^) , 

k^ 


(98) 


— i 


2r 


2 f/rn 1 C'Ttt 


where 


3e + P 2e + P 


r = 


3EL.#“’ + 2E”.1« 


6 




,( 3 , 2 ) 


6 (e -I- P) 


V 




1 C 


k^ + O {k^) , (99) 


( 100 ) 


^3 e -|- P 2 e + P^ 

These relations should be understood as generalisations of diffusion and the sound propagation in 
relativistic fluids. 

(3) 

In conformal hydrodynamics, we find that six of the A) ' coefficients give non-vanishing contribu- 

( 3 ) ( 3 ) 

tions to the linearised stress-energy tensor: from A) ^ to Ag A These transport coefficients precisely 

(■3 2"] 

correspond to the first six non-conformal u) ’ . The relevant part of the linearised conformal 
third-order stress-energy tensor (72) is then 

T+o„jin = (■'f’ + 4” + +) 8xA'8x<“«‘> + (Af + A® + Af) 8x<“8i''>8x,t.“. (101) 

(3) 

As only two combinations of the six A) enter the linearised stress-energy tensor in flat space, it 
is convenient to define the following two combinations of the six conformal transport coefficients: 


0, ^ - (Af> + Af + Af>) , 


(), = -|Af+A.f>+A, 


,(3) 


Hence, the conformal shear and sound dispersion relations become 


shear: w = -i- 


r] 


e + P 


-k^ — i 


ffTn 


1 01 


sound: oj = +Csk — iVcK =F 


(e + P) 

r, 


2c, 


2e-4P 
(Tc - 2c?Tn) k^ - i 


k^ + O {k^) , 

8 r/^rn 1 0i -|- 02 


9 (e + P)^ 3 e + P 


( 102 ) 

(103) 

(104) 

k^ + 0{k^), (105) 


where the conformal T = Tc, in the absence of bulk viscosity C; is 

r - ^ 0 

3e + P’ 

and the speed of sound c, is hxed to c, = l/\/3 for conformal fluids in four dimensions. 


(106) 
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B. Stress-energy tensor two-point function 


A particularly simple and useful two-point function of the stress-energy tensor can be 
obtained by perturbing the flat Minkowski metric tensor in the gxy component, with the momentum 
flowing in the z-direction. Perturbing to first order in 6gxy and taking a derivative with respect 
to the metric perturbation gives the scalar two-point function (T^v In this channel, the 
velocity and entropy density fields can be kept constant as they decouple from 5gxy For a detailed 
description of such a calculation, see e.g. [20, 30]. 

Using the non-conformal third-order stress-energy tensor that includes Eq. (70), we can compute 
the two-point function. We find 


T^y) = G^y’^yioj, k) = P-igu+ (gm - ^ 


i 

+ 2 


,(3,2) 

^31 


.(3,2) „,(3,2) 


-"32 


— V 


34 


UJ 


+ 2 


.(3,2) „,(3,2) 


"1 


— V. 


31 


— V. 


(3,2) 

32 




(107) 


In order to use our results for a calculation of transport coefficients in a theory with a gravitational 
dual, we need the correlator for the conformal sub-class of third-order fluids, which follows from 
using the stress-energy tensor (72). The result is 


G^y'^y^u, k) = P - iguj + [gru - 


- -Ai7 W + 


. (3) _ U3) _ U3) 
"'l "^16 


u)k^ 


(108) 


C. The Bjorken flow 

We now turn to an example of a non-linear hydrodynamic calculation in a conformally invariant 
theory: We will analyse the effects of third-order hydrodynamics on the boost-invariant Bjorken 
flow [45], which is relevant to the study of relativistic heavy ion collisions. The configuration 
describes hydrodynamic propagation of a four-dimensional boost-invariant plasma along the z- 
axis with velocity zjt. The solution can be conveniently written down in terms of proper time 
r = 'JP — the rapidity parameter ^ = arctanh(z/t) and a two-dimensional Euclidean plane, 
perpendicular to the direction of the flow. In these (flat metric) coordinates, the fluid is at rest. 
Hence, the solution for the velocity field and the metric can be written as 

= (1,0,0,0) , (109) 

gab = -dr'^ + + dx]_. (110) 

We then need to solve the hydrodynamic equations of motion for the remaining scalar field. Fol¬ 
lowing [20], it is convenient to work with the energy density e(r) instead of the entropy density 

® The name scalar is inspired by the AdS/CFT correspondence where the Sg^y ^ exp {—itui + ikz} fluctuation in the 
bulk behaves as a minimally-coupled massless scalar. Alternatively, this channel is also referred to as the tensor 
channel because the fluctuation 5gxy transforms as the spin-2 tensor fluctuation. 
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field. The function e(r) must solve Eq. (17), i.e. 


He + (e + P) = 0, (111) 

where e(r) is only a function of r. In the equation of motion, we used 11“^ as defined in Eq. (9). 
Given that the T'^ = 0 condition for a GET gives P = s/3 and that the only non-zero component 
of VaUb is = V_L^u^ = r [45], in four space-time dimensions Eq. (Ill) reduces to 

drS + ^-+TU^^ = 0. ( 112 ) 

3 T 

In conformal hydrodynamics, we have 11“^ = Using the stress-energy tensor to third order in 
the hydrodynamic expansion then gives 




At] 1 

'Sr/rn 

8 Ai' 

1 

3 r3 

9 

9 



Ap) 4Af 


2 

A3) 


+ 


3A 


(3) 


A3) 


+ 


2A 


(3) 

10 


llA 
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(3) 

11 


+ 


4A 


(3) 


5 Af^ 5A5 

+ 


(3) 


+ 


4A, 


(3) 


A3) 

'12 


6 


+ 


A3) 

U3 

6 


-A 


(3) 

15 


^ + O(r-0). (113) 


We can see that A)4 does not contribute to the energy density. Eurthermore, as discussed in Section 
(3) (3) 

III, the coefficients A)g - A20 correspond to tensors consisting solely of derivatives of the metric 
tensor that vanish in flat space. As a result, the Bjorken flow solution does not depend on any of 
those transport coefficients. 

The equation (112) can now be solved by first noting that the transport coefficients are functions 
of thermodynamical variables and thus e. Because we are working with a GET, they scale as [20] 


r] = CT] , T/TTi = Crim (^) ^ , Ai = CAi ' , A^^) ^ Qyi?) ^ (1^4) 


— — '(3) 

where C, rj, rn and A ^' are constants. The energy density function for the Bjorken flow in third- 
order hydrodynamics is then 
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A3) 

'15 
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1 

72+2u 


+ 0(r-2-3^), (115) 


where zz = 2/3. 


V. CONFORMAL TRANSPORT IN THE J\f = 4 SUPERSYMMETRIC YANG-MILLS 

THEORY 

A. Shear and sound dispersion relations 

Using our classihcation of third-order hydrodynamics, we now turn to the AdS/CET correspon¬ 
dence in order to compute the new transport coefficients 61 and 62 that entered into the linear 
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dispersion relations (104) and (105) in Sec. IV A. We will rely on the analytic results for the shear 
and the sound dispersion relations in the AA = 4 supersymmetric Yang-Mills theory at infinite 
’t Hooft coupling and infinite number of colours, Nc- The shear dispersion, computed by using 
well-known AdS/CFT techniques [25-27], was obtained in [20]. To find the dispersion of the sound 
mode, we extend the calculation of [ 20 ] to one order higher and look for the small-w and small- 
k behaviour of the lowest (hydrodynamic) quasi-normal mode of the background geometry dual 
to the AA = 4 theory at finite temperature. The geometry is the hve-dimensional, non-extremal 
asymptotically-AdS black brane solution of Einstein gravity with a negative cosmological constant: 

ds^ = [-f{u)dt^ dx^) + , (116) 

where the emblackening factor is f{u) = 1 — u^, and we have set the AdS radius to one. The 
quasi-normal mode results that we find are consistent with [40]. To quartic order in momentum k, 
the two dispersion relations are given by 


Shear: 


Sound: 


oj = — 


A-kT 
1 


k^ 


U = ±—;=k — 

Vs 


_ f (1 - hr 2 ) 4 
327r3r3 

i ^2 3 — 2 In 2 ^3 

127rr 96v^7r2r2 


i (vr^ - 24-k241n2 - 12In^ 2) 
69127r3r3 * 


(117) 

(118) 


The k'^ term of the shear dispersion relation only depends on the third-order hydrodynamic 
transport coefficient 9i, as found in Eq. (104), while the sound dispersion relation depends on 
both 9i and 92, as can be seen from Eq. (105). Hence, we can immediately determine their values 
in the M = A theory, by using the already known transport coefficients at first and second order 
[20, 21]. To date, these coefficients are known up to leading-order corrections in the ’t Hooft 
coupling, 7 = A ^/^C(3)/8, where A = gyj^Nc [36, 58-64]. The complete list was presented in [36] 
by including the value of A 2 : 


7 = [1 + 1357 + 0(7^)], 

m = 


K = 


(2-ln2) 3 ^ 

2,rT iiT ’ 

AT2rp2 

'' [ 1-107 + 0 ( 7 ^)], 


8 

l\T2rp2 

Ai = [1 + 3507 + 0(7^)] 


A2 — — 


16 
N2t2 


16 

^ 25N^T^ 2x 

As =-^ 7 + 0 ( 7 ") 


[21n2 + 5(97 + 541n2)7 + 0 ( 7 ^)] , 


(119) 

( 120 ) 
( 121 ) 
( 122 ) 

(123) 

(124) 


Using these results, along with rj/s = l/( 47 r) + O ( 7 ) [65, 66 ], and s = /2, we find the new 

third-order coefficients to be 


N'^T 

h = ^ + o(-,), 

N^T ( 7r2 


9. = 


3847r 


22 - — -181n2 + ln ^2 ) + O ( 7 ) . 


(125) 


(126) 
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B. Stress-energy tensor two-point function 

In this section, we present the scalar channel, retarded two-point function of the stress energy 
tensor G^y ;j.y{oj,k) to third order in the gradient expansion by following [25, 67]. The correlation 
function can be computed by perturbing the gxy component of the background metric tensor (116). 
The fluctuation of gxy, i.e. the gxy —>• gxy + hxy must be propagating along the z direction. We can 
then dehne (j) = uhxy/ so that (j) = h^y, giving us the momentum space equation of motion 

4>k + 4>k H-^^^2— 

where to = uj/{2ttT) and q = k/{2iTT) are the dimensionless energy-momentum variables. In order 
to recover the retarded two-point function, the solution must be in-falling at the horizon, 

0fc(tt) = (1 - Fk{u), (128) 

and Fk{u) must be regular at the horizon [u = 1). The expression for Fk{u) to second order in the 
hydrodynamic expansion, i.e. O (tu^, toq, q^), and the prescription for computing the retarded two- 
point function can be found in [25, 67]. For our purposes, we again need to extend this calculation 
to one order higher. Because the expression for Fk{u) is rather lengthy, it will not be stated here. 
We only present the two-point function with the third-order hydrodynamic corrections: 

G™(tt),q) = + (l-ln 2 )tn 2 + q2 

+ 2 In 2 - In^ 2^ - 2 In 2 iroq^ | . (129) 

Using the field theoretic result from Eq. (108), we can find two new (linear combinations of the) 
transport coefficients: 

Af ^ = ^ (^ + 42 - In^ 2) + O (7) , (130) 

C. The holographic Bjorken flow 

For our final calculation, we consider the holographic dual of the boost-invariant Bjorken flow 
in the supersymmetric W = 4 Yang-Mills theory, which was developed in [46-48, 68, 69]. The 
analytic result for the energy density written as a function of proper time, e(T), to third order in 
the gradient expansion was found in [70]. It is given by 

f 1 1 /l + 21n2\ I /3-27r2-241n2 + 241n2 2\ 1 I 

^12^3 y V 324^2 • 31/4 J ^ ^ ^ 
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Using the field theoretic result from Eq. (115) together with the knowledge of the standard N = A 
results from second-order hydrodynamics: 




TT 


■q = -N^T^ 


m = 


2-In2 
2'kT ’ 


Ai — 


V 


2-kT' 


we can first find 

r = ^ 

27r2’ 


V = 


V 2 ■ 33/4' 


Tn = 


33/^ (2-In 2) 

^71 ^ 


Ai — 


4^3’ 


(133) 


(134) 


These constants then allow us to solve for the third-order ones, giving us 


A 


(3) 


-h- 


6 + 
37^) 


4A 


(3) 




4A 


(3) 




5A 


(3) 


1(3) 




2A 


(3) 

10 


llA 


(3) 

11 


5A 


(3) 


6 




4A 


(3) 

6 


1(3) 


^12 -^13 \ (3) _ 15 - 27r^ - 45 In 2 -k 24 In'" 2 

-^ -A^5 -- 


3 

T(3) 

Ai 


6 


6 


16272 - 33 /^ 


(135) 


Using the scaling relations (114) finally gives us a linear combination of 14 third-order transport 
coefficients in the M = A theory at infinite’t Hooft coupling: 


A 


(3) 


6 




4A 


(3) 




4A 


(3) 




5A 


(3) 


q\(3) \(3) 

I _|_ ^9 


2A 


(3) 

10 


llA 


(3) 

11 


5A 


(3) 




4A, 


(3) U3) 


6 


7 


,(3) 

U2 




a(3) 

-^13 _ x(3) _ 
6 “ 


N^T 


64877 

+ 0(7). 


(15 - 277^ - 45 In 2 24 In^ 2) 

(136) 


The’t Hooft coupling corrections to all five linear combinations of the conformal third-order trans¬ 
port coefficients found in this section remain to be computed. 


VI. DISCUSSION 

In this paper, we presented a systematic algorithm for constructing the tensors that can be used 
to build the hydrodynamic gradient expansion at any order. We then used it to classify third-order 
hydrodynamics of uncharged fluids, i.e. the stress-energy tensor in the absence of Noether currents. 
The conservation equation for this stress-energy tensor therefore represents the most general, non¬ 
linear extension of the relativistic Navier-Stokes equations to the next-to-leading-order in the small 
energy-momentum gradient expansion. 

In the non-conformal case, we found that 23 scalars and 45 independent tensors could be included 
in the gradient expansion, thus giving us in total 68 new transport coefficients. In terms of the first- 
order Navier-Stokes fluids, the transport coefficients from the scalar and the tensor channels can 
be thought of as higher-order corrections to the bulk and the shear viscosity terms, respectively. 
Together with the 2 and the 15 transport coefficients in first and second-order expansions, this 
implies that there are 85 coefficients that are required to completely describe uncharged fluids to 
third order. 
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For the conformal sub-class of fluids, which have a traceless homogeneously transforming stress- 
energy tensor under the non-linear Weyl transformations, we further showed that 20 linearly in¬ 
dependent tensors could be included at third order in the gradient expansion, resulting in 20 new 
conformal transport coefficients. Hence, by including the 1 and 5 conformal transport coefficients 
in first and second-order hydrodynamics, we conclude that conformal fluids are described by 26 
transport coefficients (and independent tensor structures) to third order. 

Using our results, we computed the linear dispersion relations, uj{k), for the shear and sound 

mode, to fourth order in k. These corrections to diffusion and the sound propagation turned 

out to depend on 9 and 6 transport coefficients in the non-conformal and the conformal cases, 

respectively. By employing the AdS/CFT calculations of the dispersion relations for the AA = 4 

supersymmetric Yang-Mills theory at infinite’t Hooft coupling and number of colours, we were then 

able to determine the values of two new third-order transport coefficients in this theory, which we 

(3) 

named 6 i and 62 - Each of the two 9i depend on three of the conformal transport coefficients. A) . 
In flat space, these are the only two independent transport coefficients (or their combinations) that 
enter the linearised stress-energy tensor. Within linear response, we further used the third-order 
stress-energy tensor to compute the scalar (spin-2) two-point function, which gave us another two 
linear combinations of the transport coefficients that could be computed from holography. Lastly, 
by computing the third-order correction to the energy density of the non-linear boost-invariant 
Bjorken flow, we were able to find an additional linear combination of the new transport coefficients 
in the AA = 4 theory. Altogether, these calculations provided five linear combinations of the third- 
order transport coefficients in the AA = 4 theory at infinite ’t Hooft coupling. Another fifteen 
remain to be computed in the future; a calculation that can be done either by finding and using 
the Kubo formulae for four-point Green’s functions [30, 36, 64, 71-73] or employing the tensorial 
structure found in this paper for a calculation within the fluid/gravity correspondence [21, 74]. 

Having classified third-order hydrodynamics, numerous questions regarding its details remain 
to be answered in the future. By following the work of [30] in second-order hydrodynamics, one 
important task is to determine the relevant Kubo formulae with four-point correlation functions of 
the stress-energy tensor that would allow for a computation of all the transport coefficients, beyond 
the ones that enter into the linearised dispersion relations. These relations could then be used to 
compute the transport coefficients in perturbative and lattice QFT, as well as in strongly coupled 
held theories by using the gauge-gravity techniques. A particularly interesting task would be to 
study the huid/gravity duality [21] to third order in the gradient expansion, which could not only 
give us the microscopic values of the transport coefficients in a held theory dual to a particular 
gravitational setup, but also conhrm the counting of the linearly independent tensors identihed in 
this work. 

Perhaps the most important remaining task is the understanding of the entropy current in 
higher-order hydrodynamics [22, 23, 51, 52, 54, 75, 76]. Firstly, the precise form at second-order 
remains to be found, i.e. fully expressed in terms of the second-order transport coefficients in the 
stress-energy tensor. As argued in [22], third-order hydrodynamics may be required to resolve this 
issue even at second order. The construction of the entropy current at third order also remains to 
be performed so that potential constraints on the new maximal set of the transport coefficients can 
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be uncovered following from the non-negativity of local entropy production [22, 23, 31, 33, 51-54]. 
As a check of these relations, it would also be interesting to use the gauge-gravity duality and 
compute all non-conformal second-order transport coefficients in a field theory with a holographic 
dual. 

In gauge-gravity duality, it has been known that certain combinations of transport coefficients 
give universal values in large classes of gravitational setups. In first and second-order hydrody¬ 
namics, universality has been found for the values oi r]/s = l/(47r) and 2r]Tii — 4Ai — A 2 = 0 
[32, 36, 66, 77-79]. It is therefore natural to ask whether there exists a similar universal rela¬ 
tion between conformal (and non-conformal) third-order transport coefficients. Furthermore, as 
discussed in [36, 80], these universal relations may play an important role in the minimisation 
of entropy production in strongly coupled fluids, appearing as the coefficients of different tensor 
structures in the entropy current. 

Finally, it would be extremely interesting if some of the second and the new third-order transport 
coefficients could be measured in real fluids, either through more precise measurements of various 
dispersion relations, or by other means. In particular, the microscopic information about the fluid’s 
constituents contained in these higher order transport coefficients may play an important role in 
the modelling and studying of high-energy fluids, such as the quark-gluon plasma. 
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Appendix A: The basis of tensors for the gradient expansion 

In this Appendix, we present the proof of the claim in Section IIA that it is sufficient to use 
only transverse derivatives of and Ins, and the covariant derivatives of the Riemann tensor to 
write down the entire hydrodynamic gradient expansion at all orders. 

Let us begin by considering the action of covariant derivatives on the hydrodynamic variables, 

V l)Ua — V _i_hUa Uf, DUa , 

Vfe In s = Vj_b In s — D In s . 


(Al) 

(A2) 
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Using the equations of motion (21) and (22), we can write 

VbUa = V_L6Ua + c^UbV±a In s + higher derivatives , (A3) 

Vf,lns = V_L;,lns + V ± ■ u + higher derivatives , (A4) 

which shows that all one-derivative combinations can easily be expressed in terms of transverse 
derivatives V_L(j, at the expense of introducing higher-derivative corrections. 

Consider now a two-derivative tensor VcVx^Ua, which can be manipulated to give 


VcV J_b^o, — VJ-c^ A-b^a, '^c D\/J_b^a — ^J-c^ J-b^a, [Al, VJ_{|] Ua “b UqS/ _Lj DUa 

= Vj_cVj_bUa — c^^tcV_LfeVj_a In s — Uc [H, Vj_;,] Ua “b higher derivatives . (A5) 


The commutator appearing above can be expressed as 


[A>,V_Lfe] Ua = 


UdV^, Afee V 


Ua — D'^a '^d^be 




(A6) 


The first two terms in (A6) contain only single derivatives of Ua, so as before, the equations of 
motion allow us to eliminate the longitudinal derivatives. The commutator term can be written as 
a Riemann tensor, since for an arbitrary vector field Va, 


[Vd,Va]Vb = R%adVa. 


(A7) 


We find that 

Vo'S/±bUa = —Uc^^RfaedU^u^ + terms with (V_LcV_LbMa, V_LcV_L6 In s) + higher derivatives . 

(A8) 

Therefore, as claimed, the most general second derivative of Ua can be expressed in terms of only 
transverse derivatives, the Riemann tensor and sub-leading higher-derivative corrections, which are 
irrelevant at the order of the gradient expansion to which we are working. 

R is easy to see that with the help of the identity 

[Vrf, Vc] Vb^b2...br, = -^\icd^sfe2...6n + . .fen + • • • + Vcd^fel fe2 .. .S ) (^-9) 

this discussion can be iteratively generalised to tensors of an arbitrary order, such as 

VcV_LbiV_Lb2 .. .Ua . (AlO) 

Schematically, the proof proceed as follows: we begin by writing 

VcV• • • Ua = V_LcV±f)i^±b2 ’ ’ ’ 71VV_Lf,2 . . . Ua ■ (-^77) 

We then commute the outer two derivatives [D and Vj_ft^) in the second term, which gives us 

VcV_Lfei V_Lfe 2 . . . Ua = . . . - UcV_L6i 71V_L62 ■■■Ua- (A12) 

The ellipsis (...) in Eq. (A12) stands for various terms composed of purely transverse derivatives, 
the Riemann tensors coming from Eq. (A9) and higher-derivative terms. We must then proceed 
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by commuting D all the way through the expression (to the right), until it acts only on Ua and the 
equations of motion can be used for the last time to eliminate the remaining longitudinal derivative. 

Hence, we have shown that all of the tensorial ingredients that need to be considered in the 
construction of the hydrodynamic gradient expansion are the transverse derivatives of and Ins, 

■■■V±fe„Ua, V_Lfe^V_Lfe 2 ■■■V±fe„lns, (A13) 

their products at all possible orders, and the Riemann tensor with various metric contractions and 
arbitrary covariant derivatives acting on it, 

Vei Ve 2 • • • Ve„ Rabcd ■ (A14) 

Because transverse derivatives have the property stated in Eq. (20), the tensorial ingredients 
can be grouped into two sets: those with transverse derivatives and those with covariant derivatives, 
which always involve the Riemann tensor and its contractions. 


Appendix B: Construction of the hydrodynamics tensors that contain no Riemann tensor 

Here, we prove the claim that when constructing scalars, transverse vectors and TST two- 
tensors, we can omit from altogether and consider only the metric contractions of ±aUb 
and Vx^lns inside the hydrodynamic tensors, at all orders in the gradient expansion. 

Let us consider what happens to tensors built solely out of and V_L(jlns that get con¬ 

tracted with tt“. It is clear from Eq. (20) that for any tensor 

^^“V_LaEbi 62 ... = 0 ) (Bl) 

where Vb^b 2 ... can contain more transverse derivatives as well as combinations of and Ins. Now, 
consider a term of the form V_L(iE, 263 ...- With the help of Eq. (20), we can commute 

through the outer derivative to find 


^^“V_Lb^Vj_aEb2b3... — -g°'^'V±f,^Ua'V±cVb2b3... ■ (B2) 

The point is that the right-hand-side of Eq. (B2) could have been constructed out oiV _Lbi'^a'^±c^b 2 b 3 
in some set and gab from without any use of from By employing the same logic, 
it is easy to see that any contraction of the form 

■ ■ ■ V_LaUcEd3rf2... > (B3) 

can also be expressed as a sum of metric contractions acting on various components of No 

tensor with an undifferentiated can ever be linearly independent from tensors that use only 
transverse derivatives acting on 

By going through the same process of commuting through the expression (cf. Appendix A), 
we can also see that the same conclusion can be drawn regarding tensors of the form 

U V _Lbj V _Lb2 ■ ■ ■ UaVc\C2... ) 


(B4) 
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where now we use Eq. (41) and its derivatives, instead of (20). 

The only other possibilities that remain to be discussed are potentially uncontracted factors of 
with no derivatives acting on them. Such cases can appear in the construction of the transverse 
vectors and the TST tensors. A possible vector could have the form 

V“ = u^V, (B5) 

where all indices inside V are contracted. However, such a vector can never be transverse, i.e. 
7 ^ 0. On the other hand, a TST two-tensor with undifferentiated could have the form 

^ab^^{avb)_ (Bg) 

However, considering such tensors is redundant, as of this form could only be identically zero, 
i.e. = 0, due to various contractions of with Aab in (42). 

We have therefore shown that as claimed, only gab can be used from when constructing 
tensors that contain only transverse derivatives. As a result, hnding tensors with Rabcd^ or covariant 
derivatives Va of its various contractions, still requires us to use the full procedure outlined in Sec. 
HB. 


Appendix C: Tensors in non-conformal third-order hydrodynamics 

In this Appendix, we list the transverse, symmetric and traceless two-tensors that participate 
in the third-order gradient expansion of the non-conformal stress-energy tensor: 



'-rah 

'2 

= Vxyx‘“Vi““*> 

'-rab 

'3 

rf = 

'-rah 

'5 

= Vx‘“Vxyx'”“' 

'-rab 

>6 

Tf = Ins 

'-rah 

'8 

= Ins 

'-rab 

>9 

= V^X±"lnsV^^V> 

'-rah 

/ll 

= V/lnsV^^“v/^u^ 

'-rab 

>12 

r{§ = V/lnsV^^“V^X^ 

n^ah 

'14 

= InsV^Xy^^^"" 

'-rab 


n~'ab 

'17 

= InsV/^ InsV^X 

'-rab 

>18 

r{§ = In sV/ In 

'-Tab 

'20 

= V^^lnsV^'=lnsV^^“u^> 

'-rab 

>21 


'-rah 

'23 


'-rab 

'24 

Tg = V gu- 

'-rab 

'26 

= V gu^^V gu^'>V 

'-rab 

>27 

rg = V gu^ 

'-rah 

'29 


'-rab 

>30 

n~ab _ n,c \7 d(® 

/ 3 I - ^ c 

'-rah 

'32 

= TVaR^°-^'> 

'-rab 

>33 

Tg = 

'-rab 

'35 

= Ins 

'-rab 

'36 

Tg = 

'-rab 

'38 

= g'^'cVgu^'' 

'-rab 

>39 

rg = 

'-rab 

'41 


'-rab 

'42 

Tg = u^W^V 

'-rah 

'44 


'-rab 

'45 


Vx'°vTVxX 

Ins 

V/lnsV^^V^^V> 

InsV/V^X 

lnsV^'=lnsV^^u^> 
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Appendix D: Scalars in non-conformal third-order hydrodynamics 

The following is the list of scalars that participate in the third-order gradient expansion of the 
non-conformal stress-energy tensor: 


5i = 

52 = 

53 = 

54 = V_LX'^±bV_L'’lns 

55 = V_Lt,Vj_„lnsV/u“ 

56 = V^“lnsVj^^V/M„ 

5? = Vj^“lnsV^„V_LbM^ 

58 = V_l“ lnsV_LbVj_„u'’ 

Sg = Vj_^lnsVj_“lnsVj_f,it^ 

5io = V In sV In sVj_j'u^ 

5ii = 

5i 2 = V^„M“V_L^UbV_L‘=u'' 

5i3 = 

5i 4 = V_L^M“V_LXV_LcWb 

5i5 = 

5i6 = u-V.Rj' 

5i7 = 

5i8 = 

5i9 = jJ’lns 

520 = ±a''R 

521 = 

523 = 

522 = 


Appendix E: Vectors in non-conformal third-order hydrodynamics 

Here, we write down the list of vectors that participate in the third-order gradient expansion of 
the non-conformal stress-energy tensor: 


V“ = V/V^b'^/lns 


= V^“lnsV^bV/lns 

Vt = 

VI- 

= 

vf = 

V“ 

|/g 

= V/ix“V^^V^bix^ 

V“ = 

Pg 


Vg“ = 

*^10 

= ^±y^^±c^±bU- 

Vfi = V 

*^12 

= V^“ InsV^b InsV/ In 

V?3 = V/lnsV^XVxX 

V“ 

= V/lnsV^X'^± 6 «“ 

Vf 5 = V_|_^ In sVj_^u'^V_|_“ttb 

V“ 

= V_|_“ In sVj_^UbV_|_'^u^ 

Vfy = V_|_^ In 

V“ 

*^18 

= V_|_“ In sV_|_jtt^V_|_'^u^ 

Vfg = V_|_^ In 

V“ 

*^20 

= V_|_^ In sV_|_‘^ix“V_|_bttc 

Vfi = V_|_^ In sV_|_“ti'^V_|_bttc 

V“ 

*^22 

= 

V 23 = -RV_|_“ In s 

V“ 

*^24 

= u^u^Rfj^Vjf- In s 

Vf 5 = «Vv/lnsi?“,,, 

V“ 

•^26 


VI 7 = 

V“ 

*^28 

= u^V/u'^R^^, 
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